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First results from the Pierre Auger Observatory 
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We review in these notes the status of the construction of the Pierre Auger Observatory and present the 
first Physics results, based on the data collected during the first year and a half of operation. These results 
are preliminary, once the work to understand the systematics of the detectors are still underway. We discuss 
the cosmic ray spectrum above 3 EeV, based on the measurement done using the Surface Detector and the 
Fluorescence Detector, both, components of the observatory. We discuss, as well, the search for anisotropy near 
the Galactic Center and the limit on the photon fraction at the highest energies. 



I. INTRODUCTION 

The cosmic ray spectrum, measured at the top of the atmo- 
sphere covers a huge range in energy, going from 10 MeV to 
energies above 10^^ eV, with a differential flux that spans 31 
decades 1 1 ] . The techniques to survey this spectrum goes from 
instruments aboard satellite flights, balloon borne detectors, to 
counters that monitor the fluxes of neutrons and muons at the 
Earth surface, and at higher energies to wide area arrays of 
particle detectors. The spectrum can be divided roughly into 
four regions with very distinct behavior. The first one, with 
energies below 1 GeV, has a very distinctive character from 
the rest of it. Its shape and cut-off is strongly dependent on 
the phase of the solar cycle, a phenomenon known as solar 
modulation. Actually, there is an inverse correlation between 
the intensity of cosmic rays at the top of the atmosphere and 
the level of the solar activity 

The region above 1 GeV show a spectrum with a power 
law dependence, N{E) dE = K dE, where the spectral 
index x varies as 2.7 < x < 3.2. The region between 1 GeV 
and the knee region at 4 x 10^^ eV, is characterized by an in- 
dex X :^ 2.7. These cosmic rays most likely are produced 
at supernova explosions and their remnants |4]. At the knee 
(4 X 10^^ eV) the power law index steepens to 3.2 until the so 
called ankle, at 5 x 10^^ eV. The origin of the cosmic rays in 
this region is less clear and subject of much conjecture. Above 
the ankle the spectrum flattens again to an index x :^ 2.8 and 
this is interpreted by many authors as a cross over from the 
steeper galactic component to a harder extra galactic source 
for the cosmic rays QEh- 

The existence of cosmic rays with energies above 10^^ eV 
presents a puzzling problem in high energy astrophysics. The 
first event of this class of phenomena was observed at the be- 
ginning of the 60's by John Linsley at the Volcano Ranch ex- 
periment, in New Mexico Since then, many events were 
detected, at different sites, using quite distinct techniques 
[ni uA ' If those cosmic rays are common mat- 
ter, that is, protons, nuclei or even photons, they undergo well 
known nuclear and electromagnetic processes, during their 
propagation through space. Their energies are degraded by 
the interaction with the cosmic radiation background (CMB), 
way before they reached the Earth | ij, [l^. After travelling 
distances at the scale of 50 Mpc, their energies should be un- 



der 10^^ eV, thus restricting the possible conventional astro- 
physical objects, which could be sources of them, and those 
should be easily localized through astronomical instruments. 
On the other hand, it is very difficult to explain the acceler- 
ation of charged particles, with energies up to 3 x 10^^ eV 
or even greater, on known astrophysical objects, by means of 
electromagnetic forces, the only conventional ones capable of 
long range and long periods of acceleration. 



II. THE PIERRE AUGER OBSERVATORY 

The Pierre Auger Observatory (3, [3 was designed to 
study the higher - above 10^^ eV - end of the cosmic ray 
spectrum, with high statistics over the whole sky. The detec- 
tors are optimized to measure the energy spectrum, the direc- 
tions of arrival and the chemical composition of the cosmic 
rays, using two complementary techniques, surface detectors 
(SD) based on Cherenkov radiators and fluorescence light de- 
tectors (ED). The complete project calls for two sites, one in 
the southern hemisphere, well ahead in its construction, and 
a northern site, in order to achive the homogeneous coverage 
of the whole sky, essential to pinpoint the sources of the ultra 
high-energy cosmic rays. 

The southern site of the Observatory, located in the 
province of Mendoza, in Argentina, at the latitude 35° South 
and longitude 69° West, on a very flat plateau at 1 400 m above 
sea level, is bounded by the Andes Cordillera on the west. The 
region has a very clear sky, with little light pollution, essential 
for the operation of the fluorescence system. The main office 
of the observatory is located at the northern entrance of the 
City of Malargiie. The construction of this part of the detector 
should be completed on the second half of 2007. The layout of 
the site is shown in Eigure[ll indicating the distribution of de- 
tectors already deployed and the position of the fluorescence 
buildings. 

The construction of the northern site, in the southern part of 
the state of Colorado, at the United States, will start after the 
full completion of the southern site. 

The initial phase of the Pierre Auger Observatory, the Engi- 
neering Array was built during the period 2000-2001. Then, 
40 Cherenkov tanks were laid and instrumented, to test the 
components of the array and to prove the soundness of its de- 
sign. They were set in a roughly hexagonal array, covering an 
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FIG. 1 : Layout of the Pierre Auger Observatory southern site show- 
ing the four FD telescopes (eyes). The dots represents the deploye- 
ment location of the SD stations. The darker shade represents the 
area instrumented in July 2006. The FD eyes Los Leones, Los Mora- 
dos and Coihueco are in full operation, while the construction of 
Loma Amarilla and its instrumentation should be completed by the 
end of 2006. 



on flatbed trucks, carrying six tanks on a load. After prepa- 
ration which includes the mounting of the electronic compo- 
nents, the cabling and fitting the internal Tyvek bag, the tanks 
are deployed on the field, on prepared ground, using 4-wheel 
drive vehicles. 
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area of 54 km^. Two of the tanks, at the center of the array, 
were laid side by side, in order to cross calibrate their signals. 
Two prototype telescopes were installed on the Los Leones 
hill, overlooking this ground array. They were successfully 
tested during the southern summer of 2001-2002. The lessons 
learned in the Engineering Array lead to improvement in the 
final design of the detector components il8>1 . 



A. The Surface Detector 

The Surface Detector (SD) is a ground array which spans 
an area of 3 000 km^, with 1 600 stations set on a regular 
triangular grid, with 1 500 m separation between them. The 
stations communicate with the central base station through a 
radio link. 

Each station at the SD is a cylindrical tank, filled with 
12 000 1 of purified water, operating as a Cherenkov light 
detector. The tanks are manufactured by rotational molding 
process, using a high density polyethylene resin, with 12.7 
mm thick walls, opaque to external light. The walls have two 
layers, an external in the color beige, to minimize the environ- 
mental impact and reduce the heat absorption, and the internal 
black, with the addition of carbon black to the resin. The wa- 
ter is contained in a liner inside the tank, a bag made of a 
sandwich of polyolefin-Ty vek film. The Tyvek film has a high 
reflectivity to ultraviolet light and its role is to diffuse the UV 
Cherenkov light within the volume of water. 

The tanks are prepared in the Assembly Building, just be- 
hid the Main campus building. The water purification plant 
produces enough water for 3 tanks per day, with a resistivity 
of 15 MQ, per cm. The tanks are manufactured in Sao Paulo, 
Brazil, and Buenos Aires and shipped all the way to Malargiie 



FIG. 2: A schematic view of the Cherenkov water tanks, with the 
components indicated in the figure. 

The Cherenkov light, diffused within the volume of water, 
is collected by three 20.3 cm diameter photomultiplier tubes 
(Photonis XP1805), set in a symmetric pattern on top of the 
tanks, facing downwards. This arrangement avoids the direct 
hit of the Cherenkov light, collecting a signal which is ho- 
mogeneous and proportional to the length of charged tracks 
crossing the water. A schematic view of a SD station is shown 
in Figure 121 

The electronics for the detectors are housed in a box out- 
side the tanks and communicate with a base station through a 
radio WLAN, operating in the 915 MHz band. The stations are 
powered by a bank of two special 12 V batteries, which are 
fed by two solar panels of 55 W each. The time synchroniza- 
tion of the tanks is based on a GPS system, capable of a time 
alignment precision of about 10 ns 1 19]. A 7 GHz microwave 
backbone links the base stations to the central data acquisition 
station (CDAS), at the campus of Malargiie in the southern 
site. 

Each detector station has a two level trigger, a hardware im- 
plemented Tl and a software T2. The Tl trigger is decided in 
a PLD (Programmable Logical Device) 1 20], set with a thresh- 
old defined in terms of a vertical equivalent muon (VEM) 
crossing a tank, with a typical value of 1.75 VEM on a single 
25 ns time bin, in coincidence at all working PMT's. This trig- 
ger has a rate of about 100 Hz and is necessary to detect fast 
signals associated to the muons of horizontal showers. There 
is a second condition for Tl, which is a time over threshold 
condition (ToT), requiring that the signal in 13 FADC bins 
out of a window of 120 bins are above a value of 0.2 VEM. 
The ToT trigger rate averages 1.6 Hz, dominated by coinci- 
dences of two muons and is efficient to select small signals 
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away from the core of showers at the tail of the lateral dis- 
tribution. Tl is adjusted so that its rate is about 100 Hz. T2 
limits the trigger rate at each station to less than 20 Hz, so 
as not to saturate the radio bandwidth available. ToT triggers 
are promptly promoted to T2, while the threshold Tl are re- 
quested a tighter condition of 3.2 VEM in coincidence for the 
3 PMTs. The event trigger (T3) is set at CDAS combining the 
triggers of contiguous individual stations. The requirement in 
the number of stations triggered sets the lower energy thresh- 
old. Typically, four stations are required for a threshold of 
10^^ eV. The communication between the central campus and 
the individual stations is bi-directional to allow the T3 trigger 
request the data to be downloaded to the CDAS system. 
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FIG. 3: Typical histogram used in the calibration of the surface sta- 
tions. This was taken from a prototype tank equipped with scintilla- 
tion pads to select the component of the signal due to vertical muons. 
The hump of the experimental distribution of single muons is related 
to the signal due to a vertical incoming muon 1 21]. 



The calibration of the surface stations is done continuously. 
At every minute the histogram of low energy particles is taken, 
correspondinding roughly to about 100 000 events, mostly at- 
mospheric muons coming from all angles into the detector. 
This calibration is performed in parallel to the data taking and 
at every 6 minutes this is sent to the central data acquisision 
system (CDAS) for monitoring |21]. Figure |3l shows a the 
histogram for a surface station, taken at a prototype detector. 

The incoming angle of a shower is reconstructed from the 
timing of arrival of the signals in the tanks. The detectors have 
a large cross section even for a shower with very hi gh z enith 
angle, so that Auger is quite sensitive to neutrinos I22l 12311 . 
A highly inclined shower, originating from a hadronic cosmic 
ray, has a very characteristic signature, having lost a substan- 
tial part of its electromagnetic component, with only a core of 
energetic muons remaining. The muons arrive concentrated 
in time, generating signals with a very sharp peak. In con- 
trast, a shower which has a large electromagnetic component 
is spread in time, with a much more complex structure. 



B. The Fluorescence Detector 

The Fluorescence Detector (FD) is composed of 4 eyes dis- 
posed on the vertices of a diamond-like configuration, with 
a separation of 65.7 km along the axis South-North (actually 
this axis is tilted by 19° into Northeast) and 57.0 km in the 
West-East direction (tilted by 20° into Southeast). They are 
at the periphery of the ground array and all stations of the SD 
are contained in the field of view (FOV) of the FD telescopes 
1 24] (see FigurelH). 




FIG. 4: The Coihueco telescope building, with the shutters closed. 
Lower picture: The diaphragm of a telescope, with corrector rings. 
The picture shows the camera, with its 440 pixels facing the mirror. 

Each eye has six independent telescopes, each with a field 
of view of 30° in azimuth and 28.6° in elevation, adding to 
a 180° view of the array. Figure |4| shows the building of the 
Coihueco eye located on the west side of the array (see Fig- 
ure [TJ. The fluorescence light is collected by a mirror with 
a radius of 3.4 m and reflected into a camera, located at the 
focal surface of the mirror. The telescopes use a Schmidt op- 
tics design to avoid coma aberration, with a diaphragm, at the 
center of curvature of the mirror, with an external radius of 
0.85 m. The advantage of this scheme is to project a point 
of light in the sky into a reasonably homogeneous spot, at the 
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focal surface, of 0.25° radius. 

The shape of the telescope mirror is a square with rounded 
corners with a side of 3.8 m. The radius of the diaphragm may 
be enlarged to 1 . 1 m, by adding a corrector lens annulus on the 
area with radius between 0.85 m and 1.1 m |25] (see Figure 
|4|). The telescope, which is remotely operated, is protected 
from the environment by an external shutter and a MUG- 6 
UV transmitting filter, which reduces the light background. 

The camera is composed by 440 pixels, each a hexagonal 
photomultiplier (Photonis XP3062), which monitors a solid 
angle of (1.5°)^ projected onto the sky. The dead spaces be- 
tween the photomultipliers is corrected by a reflective wedge 
device, called the Mercedes corrector, making the sky expo- 
sure very uniform. 

The electronics of the FD detector was designed to be op- 
erated remotely and with flexibility to reprogram the trigger 
to accommodate non-standard physical processes which may 
show up. The trigger rate is not limited by the transmission 
band available. The FD communicates with the CDAS sys- 
tem through the 7 GHz microwave backbone, with a capacity 
of 34 Mbps throughput. 

To be able to operate in hybrid mode the absolute time 
alignment with the Cherenkov water detector must be better 
than 120 ns. The trigger for each FD pixel is programmed to 
have a rate below 100 Hz, while a special processor, with a 
built-in pattern recognition algorithm, keeps the overall trig- 
ger rate below 0.1 Hz. The pixels are sampled and the signal 
digitized at every 100 ns. The first level trigger uses a boxcar 
type addition of the signal over a 1 /iS bin and a pixel trigger 
set whenever this sum is 3 a above the background noise. An 
event is triggered whenever a set of 5 contiguous pixels are 
triggered with a time sequence associated |26]. 

The calibration of the camera is performed by exposing the 
PMTs to signals from calibrated light sources diffused from an 
aparatus mounted on the external window of the telescopes. 
The camera is illuminated uniformly and the gain of each 
PMT characterized OEjl. 

The determination of the shower energy requires an ac- 
curate estimate of the atmosphere attenuation, due to the 
Rayleigh (molecular) and Mie (aerosol) scattering, of the light 
emitted by the fluorescence excitation of the atmosphere. This 
is done by a complex set of detectors, the Horizontal Attenu- 
ation Monitor (HAM), the Aerosol Phase Function monitors 
(APF) and Lidar systems mounted at each eye. The 

system is complemented by sky monitoring CCD's, which 
measures the light output of selected stars. A set of infrared 
cameras at each eye monitors the cloud coverage at the site. 
Monitoring the background noise due to stars in the field of 
view of the pixels are tools used to get the proper alignment 
of the telescopes. 

The Central Laser Facility (CLF) 1 3 1 ] is a steerable auto- 
matic system which produces regular pulses of linearly polar- 
ized UV light at 355nm. It is located in the middle of the array, 
26 km away from Los Leones and 34 km from the Coiheco 
eye. This system provides a complementary measurement of 
the aerosol vertical optical depth versus height, and the hori- 
zontal uniformity of the atmosphere across the aperture of the 
array. This system creates an artificial cosmic ray by feeding 



a signal into a nearby tank through a fiber optics cable. The 
laser track in the atmosphere is intense enough to be registred 
by all fluorescence detectors. The time recorded at each de- 
tector can be used to measure and monitor the relative timing 
between the SD stations and the FD eyes. This time offset has 
been measured to better than 50 ns |32], leading to a system- 
atic uncertainty in the core location of 20 m. This is shown in 
Figure 13 
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FIG. 5: Resolution of the reconstructed position for vertical laser 
shots from the CLF. 

Due to the much improved angular accuracy, the hybrid 
data sample is ideal for anisotropy studies. Many ground pa- 
rameters, like the shower front curvature and thickness, have 
always been difficult to measure experimentally, and were 
usually determined from Monte Carlos simulation. The hy- 
brid sample provides a unique opportunity in this respect. As 
mentioned, the geometrical reconstruction can be done using 
only one ground station, thus all the remaining detectors can 
be used to measure the shower characteristics. The possibility 
of studying the same set of air showers with two independent 
methods is valuable in understanding the strengths and limita- 
tions of each technique. The hybrid analysis benefits from the 
calorimetry of the fluorescence technique and the uniformity 
of the surface detector aperture |33]. 

We show in the Figure|6|a typical particle distribution at the 
ground, for a shower with an energy of 10 EeV, with a zenith 
angle of 45°. Although gammas and electrons and positrons 
dominate the shower, at large distances from the core and at 
large angles, muons come to dominate the particle distribu- 
tions. 



III. AUGER DATA 

Since the begining of operations in 2001, Auger has been 
taking data continuously, first in the Engineering Array mode, 
which is already reported 1 18]. Since January 2004 it is tak- 
ing data in Physics Mode, with results that will be reported 
in the following. A typical example of an event recorded by 
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FIG. 6: Particle distribution at ground from a shower of 10 EeV at an 
angle of 45o. 



the SD system is shown in the Figure The energy of the 
shower is infered from the lateral distribution function (LDP), 
fitted from density of particles hitting each tank at distinct dis- 
tances The density of particles at 1 km from the core, 
S(IOOO), is quite independent of the nature of the primary cos- 
mic ray, according to distributions extracted from shower sim- 
ulation programs, as mentioned before. Each station in the 
array is continuously calibrated by collecting at regular inter- 
vals signals generated by single muons. From this signal the 
value of a VEM (vertical equivalent muon) is extracted, tak- 
ing into account variations in the atmospheric temperature and 
pressure. 

[ Ljl*rmf diBlrlbutlon funclloo in | 
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FIG. 7: A typical event registered by the SD system, shown in the 
Event Display. The shower has a zenith angle of 30° and 10 stations 
are used in the fitting of the lateral distribution function. The signal 
is quoted in VEM (Vertical Equivalent Muon). 

The Auger ground stations are sensitive to very inclined 
showers once it has a reasonable cross-section due to the 1.2 m 
column of water. The very inclined showers will have trans- 
versed a larger amount of atmospheric matter before hitting 
the station and a large part of the electromagnetic shower will 
have dissipated, allowing for a much larger relative muon 
component. An example of this class of event is displyed 



in Figure [8] with a shower hitting 31 stations, coming with 
a zenith inclination of 8 8 ° . 
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FIG. 8: Event with a large zenith angle, 88° in this case. Most of 
the signal is carried by the surviving muons, once a large part of the 
electromagnetic shower has been dissipated. 

The structure of the FD events is exhibited in Figure |9l 
taken from the event monitoring system. On both pictures 
the map of the pixels are shown as if one would be facing the 
sky. So, the first event, a typical low energy shower, is flying 
top-down, right to left. The histogram on the top-right of the 
display shows the time structure of the signal in a selection of 
the pixels. The colour code is there just to correlate the pixel 
and the signal. Each bin on the right plot corresponds to 100 
ns, so that the span of the signal from 260 to 360 units, actu- 
ally corresponds to a time interval of 10 jus. The bottom picture 
of Figure|9lrepresents a typical event propagating downwards, 
while the top part there is an laser shot event, with a time struc- 
ture going upwards. 

To reconstruct the geometry of the shower, first the shower 
detector plane (SDP) is infered by optimizing the line of light 
crossing the camera, using the signals as weight. The best 
estimate of the normal vector to the SDP, fisDP, is obtained by 
minimizing 



where the signal measured in pixel / is used with the weight Wi 
and the n corresponds to the direction pointing to the source 
in the sky. The three dimensional geometry is recovered using 
the angular velocity of the signal. For each shower pixel / the 
average time of the arrival of the light at that pixel field of 
view, ti, is determined from the FADC traces. The expression 
13^, 



Rp 

U = ^0 + — tan 
c 



(Xo-li) 



allow for a fitting of the shower parameters, Rp, Xo and ^o- 
Here, c is the velocity of light, Rp the shower distance of 
closest approach to the detector and the time at which the 
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FIG. 9: Typical FD events. The camera display represent the shower 
projected on the sky. The color code just identify the pixel. The black 
dots on the pixels mark those exhibited on the right. The time axis on 
the right side histogram is set in units of 100 ns, the label 350 means 
35 ^s. Top: Typical Laser shot event, where the light goes from the 
ground up. Lower: Regular event. 

shower point reaches the position of closest approach. in- 
dicated on Figure [10| is the direction of the pixel / projected 
onto the S DP and Xo is the angle between the shower axis and 
the direction from the detector to the shower landing point. 
Figure [H] shows an example of time versus angle plot for a 
shower seen in stereo mode. However, this procedure is not 
free of ambiguities, which can be resolved with the input from 
the SD system. The timing information and location from a 
station closest to the shower landing point can be related to 
the time to, 

^tank ' ^shw 

to — Hank 7 



where Rtank is the vector connecting the fluorescence detector 
to the ground station and Sshw is the unit vector associated to 
the shower propagation. 




FIG. 10: Shower detector plane. 



To measure the energy, the light emitted by the source is 
reconstructed making the corrections for the atmosphere at- 



tenuation and than subtracting the Cherenkov component of 
the signal, identifying the fluorescence component. The time 
profile and the longitudinal profile for a stereo event event is 
exhibited in Figure [H] for both views of the event. The line 
fitting the longitudinal profile represents the Gaisser-Hillas 
function llsvll . 




FIG. 1 1 : Event seen in stereo mode. The left hand side is the view 
from Los Leones and the right from Coihueco. The lower plots shows 
the time versus angle correlation as seen by each eye. The fit pa- 
rameters are shown in the plots. Lower: The time profile and the 
longitudinal profile of the event above 



Although there is much work to do to improve the quality 
of the measurements done with the SD and FD components, 
there are some preliminary tests that points to overall quality 
of the data. In particular, the correlation in the estimation of 
the energy of hybrid events, measured by the SD and the FD, 
is quite consistent. 
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IV. AUGER RESULTS 
A. A first estimate of the cosmic ray spectrum above 3 EeV 

For the first attempt to estimate the cosmic ray spectrum 
above 3 EeV the analysis was based on a selection of events 
collected from January 1st, 2004, through June 5th, 2005. 
The event acceptance criteria and exposure calculation are de- 
scribed in a separate paper Events included in the sample 
have zenith angles of less than 60° and energy greater than 3 
EeV, resulting in a selection of 3525 events. The array is fully 
efficient for detecting such showers, so the acceptance at any 
time is given by the simple geometric aperture. The cumu- 
lative exposure adds up to 1750 km^ sr yr, which is about 
7% greater than the total exposure obtained by AGASA |39], 
along the years. The average array size during the time of this 
exposure was 22% of what will be available when the southern 
site of the Observatory has been completed. 

Assigning energies to the SD event set is a two-step pro- 
cess, of which the first is to assign an energy parameter S^^ 
to each event. Then the hybrid events are used to establish 
the rule for converting 838 to energy. The energy parameter 
^38 for each shower comes from its experimentally measured 
5'(1000), which is the time-integrated water Cherenkov signal 
that would be measured by a tank 1000 meters from the core. 

The signal 5(1000) is attenuated at large slant depths. Its 
dependence on zenith angle is derived empirically by exploit- 
ing the nearly isotropic intensity of cosmic rays. By fixing 
a specific intensity /q (counts per unit of sin^O), one finds 
for each zenith angle the value of 5'(1000) such that /(> 
5'(1000)) = Iq. We calculated a particular constant intensity 
cut curve CIC(0) relative to the value at the median zenith an- 
gle (0 ^ 38°). Given 5(1000) and for any measured shower, 
the energy parameter Sss is defined by = 5'( 1000) / C/C(0) . 
It may be regarded as the 5(1000) measurement the shower 
would have produced if it had arrived 38° from the zenith. 

The ^38 parameter is well correlated to the energy extracted 
from well reconstructed hybrid showers, as can be seen at Fig- 
ure [T2| The fitted line gives an empirical rule for assigning 
energies in EeV, based on the ^38 parameter, given in VEM, 



E = 0.16 X Sli^ = 0.16 X [5(1OOO)/C/C(0)] 



1.06 



On the top picture of Figure [2] we show the preliminary 
spectrum extracted from the data up to June 2005. Plotted on 
the vertical axis is the differential intensity 

dl dl 

=E — . 

dlnE dE 

Error bars on the points indicate the statistical uncertainty (or 
95% CL upper limit). Systematic uncertainty is indicated by 
double arrows at two different energies. On the bottom of 
Figure El the percentage deviation from the best-fit power 
law, 100 X {{dI/d{lnE) - F) /F, is shown. The fitted function 
is 



F = 30.9 ± 1.7 X (E/EeV)' 



-1.84±0.03 
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FIG. 12: Shower energy for hybrid events as a function of the surface 
detector parameter S^^. 
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FIG. 13: Top: Estimated spectrum. Lower: Deviation of the be- 
haviour of the spectrum, in relation to a power law (see text). 



with a chisquare per degree of freedom of 2.4. 
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V. ANISOTROPY STUDIES AROUND THE GALACTIC 
CENTER 

The AGASA experiment EEl and the SUGAR experi- 
ment 1 42] have reported an excess of cosmic rays coming from 
the general direction of the center of our galaxy, with energies 
in the EeV range. There is a very massive black hole at the 
center of the galaxy (CG), which provides a natural candidate 
for a cosmic ray accelerator to very high energies. Thus this 
region provides an attractive target for anisotropy studies with 
the Pierre Auger Observatory. 

For this analysis a data set was selected comprising all the 
SD events collected from January 1st, 2004, until June 6th, 
2005, with a zenith angle less than 60°. Events from the Sur- 
face Detector that passed the 3 -fold or the 4-fold data acquisi- 
tion triggers and satisfyed our high level physics trigger (T4) 
and our quality trigger (T5) |20] were selected. The T5 selec- 
tion is independent of energy and ensures a better quality for 
the event reconstruction. This data set has an angular resolu- 
tion better than 2.2° l43ll for all of the 3-fold events, regardless 
of the zenith angle considered, and better than 1.7° for events 
with multiplicities larger than 3 SD stations. 




FIG. 14: Lambert projections of the galactic centre region, GC 
(cross), galactic plane (solid line), regions of excess of AGASA and 
SUGAR (circles), AGASA f.o.v. limit (dashed line). The explana- 
tion of each figure is detailed in the text. 

To estimate the coverage map, needed to construct excess 
probability maps, a shuffling technique was used. Figure fT4K 
shows the coverage map obtained from the SD sample in a 
region around the GC. It has the same color scale as the sig- 
nificance maps, but limited to the range to 1. Figures fT4^. 
C and D, present the chance probability distributions |44] - 
mapped to positive Gaussian significance for excesses and 
negative for deficits - in the same region for various filter- 
ing and energy cuts corresponding to the various searches. 
Figures fT4b is the significance map in the range from 0.8 to 
3.2 EeV, smoothed using the individual pointing resolution of 
the events and a 1.5° filter (Auger like excess), while fT4lZ! 
is the same, smoothed at 3.7° (SUGAR like excess). Figure 
fT4b covers the range from 1.0 to 2.5 EeV, smoothed at 13.3° 
(AGASA like excess). In these maps the chance probability 
distributions are consistent with those expected as a result of 



statistical fluctuations from an isotropic sky. 

In the region where the AGASA excess was reported, there 
are 1155 events in the energy range from 1.0 to 2.5 EeV 
observed at the Auger Observatory, which is to be com- 
pared to an expected sample of 1160.7 events, thus a ratio of 
1.00 ±0.03. 

These results are not compatible with the excess reported 
by AGASA. The level of 22% excess that they observed would 
translate into a 7.5 a excess. At the worst case scenario, where 
the source would be made of nucleons and the bulk of the 
background at this energy range made of much heavier nuclei 
(e.g. Iron), the difference in detection efficiency of the Auger 
trigger at 1 EeV would reduce the sensitivity to a source ex- 
cess. However, using the ratio of Fe (70 %) to proton (50 %) 
efficiency at 1 EeV (1.44, an upper bound in the range from 1 
to 2.5 EeV) a 5.2 a event excess would still be expected from 
our data set. 

We observe in the angular/energy window, where SUGAR 
claims an excess, 144 events observed against 150.9 expected, 
with a ratio 0.95 ± 0.08. With over an order of magnitude 
more statistics we are not able to confirm this claim. 

A search was performed for signals of a point-like source 
in the direction of the galactic center using a 1.5° Gaussian 
filter corresponding to the angular resolution of the SD |43]. 
In the energy range of 0.8 to 3.2 EeV, we obtain 24.3 events 
observed, against 23.9 expected, with a ratio 1.0 ±0.1. A 
95% CL upper bound on the number of events coming from 
a point source in that window is = 6.7. This bound can be 
translated into a flux upper limit (<l>^) integrated in this energy 
range. In the simplest case in which the source has a spectrum 
similar to the one of the overall CR spectrum (&N /&E oc E~^), 

where a is the size of the Gaussian filter used. Using 

^cr{E) = 1.5 ^ {E/EeVy^ x 10"^^ (EeV-^m-^s-^sr-^), 

where ^, a number in the range from 1 to 2.5, denotes our 
uncertainty on the CR flux (^ is around unity for Auger and 
2.5 for AGASA), introducing the factor £ for the iron/proton 
detection efficiency ratio (1 < 8 < 1.6 for E in the range from 
0.8 to 3.2 EeV) and, integrating in that energy range we obtain 

<J>, < 2.6 ^ 8 X lO-^^-V^ (95% CL). 

In the worst case scenario, where both ^ and 8 take their 
maximum value, the bound is = 10-6 x 10~^^ m~^ s~^, 
excluding the neutron source scenario, suggested by the ref- 
erences El to account for the AGASA excess, or by 
in connection with the HESS measurements. More 
details about the GC anisotropy studies with the Auger Obser- 
vatory data can be found in 14811 . 



VI. AN UPPER LIMIT ON THE PRIMARY PHOTON 
FRACTION 

The photon upper limit derived here is based on the direct 
observation of the longitudinal air shower profile and makes 



Brazilian Journal of Physics, vol. 36, no. 4A, December , 2006 



9 



use of the hybrid detection technique and Xmax is used as dis- 
criminant observable. The information from triggered sur- 
face detectors in hybrid events considerably reduces the un- 
certainty in shower track geometry. 

The data were taken by the 12 fluorescence telescopes at 
the Los Leones and Coihueco eyes, from January 2004 up 
to April 2005. The number of deployed surface detector sta- 
tions grew from about 200 to 800 during this period. For the 
analysis, hybrid events were selected, i.e. showers observed 
both, by the surface tanks (at least one) and the telescopes 
^33]. Even for one triggered tank only, the additional timing 
constraint allows a significantly improved geometry fit to the 
observed profile, which leads to a reduced uncertainty in the 
reconstructed Xmax- 

The reconstruction is based on an end-to-end calibration of 
the fluorescence telescopes |281, on monitoring data of local 
atmospheric conditions |30, 491, and includes an improved 
subtraction of Cherenkov light |50] and reconstruction of en- 
ergy deposit profiles for deriving the primary energy. In total, 
16 events with energies above 10^^ eV are selected. The total 
uncertainty ^^ax ^^e reconstructed depth of shower max- 
imum is composed of several contributions which, in general, 
vary from event to event. A conservative estimate of the cur- 
rent Xmax uncertainties gives AXj^^^ :^ 40 g cm~^. Among the 
main contributions, each one in general well below AXmax = 
15 g cm~^, are the statistical uncertainty from the profile fit, 
the uncertainty in shower geometry, the uncertainty in atmo- 
spheric conditions such as the air density profile, and the un- 
certainty in the reconstructed primary energy, which is taken 
as input for the primary photon simulation. For each event, 
high-statistics shower simulations are performed for photons 
for the specific event conditions. A simulation study of the de- 
tector acceptance to photons and nuclear primaries has been 
conducted. For the chosen cuts, the ratio of the acceptance 
to photon-induced showers to that of nuclear primaries (pro- 
ton or iron nuclei) is £ = 0.88. A corresponding correction 
is applied to the derived photon limit. Figure [21 shows as an 
example an event of 16 EeV primary energy observed with 
^max = 780 g cm~^, compared to the corresponding Xmax 
distribution expected for simulated primary photons. With 
< Xmax > = 1020 g cm~^, photon showers are on average ex- 
pected to reach maximum at depths considerably greater than 
observed. Shower- to- shower fluctuations are large due to the 
LPM effect (rms of 80 g cm~^) and well in excess of the mea- 
surement uncertainty. For all 16 events, the observed Xmax is 
well below the average value expected for photons. The Xmax 
distribution of the data is also displayed in Figure [21 More 
details about this analysis can be found in 1 5 1 ] . 

The statistical method for deriving an upper limit follows 
that introduced in |52]. For the Auger data sample, an upper 
limit on the photon fraction of 26% at a confidence level of 
95% is derived. In Figure[T2l this upper limit (Auger) is plot- 
ted together with previous experimental limits, from AG AS A 
(Al) 153] , (A2) 152] and Haverah Park (HP) |54] data and 
compared to some estimates based on non-acceleration mod- 
els (ZB, SHDM and TD from |55] and SHDM' from |56]). 
The presented 26% limit confirms and improves the exist- 
ing limits above 10^^ eV. More details of this analysis can 



be found in fsT] and references therein. 
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FIG. 15: Example of Xmax measured in an individual shower of 16 
EeV (point with error bar) compared to the Xmax distribution ex- 
pected for photon showers (solid line). Also shown the Xmax dis- 
tribution of the data sample (dashed line). 
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FIG. 16: Upper limits on cosmic-ray photon fraction compared with 
previous experiments and estimates based on non- acceleration mod- 
els (see main text). 



VII. SUMMARY 

The Pierre Auger Observatory is still under construction 
but has already the largest integrated exposure to high en- 
ergy cosmic rays. The combination of fluorescence and the 
surface detector measurements allow for the reconstruction of 
the shower geometry and its energy with much greater qual- 
ity than what could be achieved with either detector standing 
alone. Each of the detectors have different systematics, al- 
lowing for valuable information for cross-checking the results 
from each of them. 

The observatory should be finished by mid 2007, accumu- 
lating by then a much larger exposure than what was used for 
the preliminary results presented here. That will allow for the 
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search of anisotropies in the southern sky, as well as the test of the predicted GZK suppression. 
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